Adiabatic Splitting, Transport, and Self- Trapping 
of a Bose-Einstein Condensate in a Double- Well Potential 
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We show that the adiabatic dynamics of a Bose-Einstein condensate (BEC) in a double well po- 
tential can be described in terms of a dark variable resulting from the combination of the population 
imbalance and the spatial atomic coherence between the two wells. By means of this dark variable, 
we extend, to the non-linear matter wave case, the recent proposal by Vitanov and Shore [Phys. 
Rev. A 73, 053402 (2006)] on adiabatic passage techniques to coherently control the population 
of two internal levels of an atom/molecule. We investigate the conditions to adiabatically split 
or transport a BEC as well as to prepare an adiabatic self trapping state by the optimal delayed 
temporal variation of the tunneling rate via either the energy bias between the two wells or the 
BEC non-linearity. The emergence of non-linear eigenstates and unstable stationary solutions of 
the system as well as their role in the breaking down of the adiabatic dynamics is investigated in 
detail. 

PACS numbers: 03.75.Kk, 03.75.Lm 



I. INTRODUCTION 

Bose Einstein condensates (BEC) in double well poten- 
tials have drawn a lot of attention both theoretically and 
experimentally for the possibilities they offer to study 
fundamental quantum mechanical effects at the macro- 
scopic level as well as for potential applications like inter- 
ferometry, high precision measurements or thermometry 
Most of the experimental realizations lead to the 
splitting of the condensate into two independent parts. 
Nevertheless, recently, also weakly linked parts of a BEC 
in a double well potential forming a single Joscphson 
junction have been achieved Q. In contrast to Joseph- 
son junctions realized in superconductors and supcrfluids 
in BEC the non-linear interatomic interactions play 
a crucial role. In the presence of the non-linearity two 
dynamical regimes have been predicted: (i) anharmonic 
Josephson oscillations Q , if the initial population imbal- 
ance of the two wells is below a critical value and (ii) 
macroscopic quantum self-trapping i.e, inhibition of 
large amplitude Joscphson oscillations above the thresh- 
old for the population imbalance. This threshold corre- 
sponds to the population imbalance for which the dif- 
ference between the two on-site interaction energies be- 
comes larger than the tunneling energy splitting. Both 
dynamical regimes have been explored experimentally in 
a single Joscphson junction and in arrays 

Recently, a lot of attention has been devoted to explore 
techniques to coherently control the non-linear dynamics 
of a BEC in double well potentials by modulating in time 
either the potential Q or the non-linearity The two 
main studied regimes are the non-linear Landau-Zener 
and the Rosen-Zener regimes. In the former 
the tunneling rate is fixed while the energy bias is linearly 
varied in time, in the later the energy bias is varied in 
time while the tunneling rate is switching on/off following 
e.g., a temporal Gaussian profile. Both these regimes 



have been deeply investigated in double [l0l - [T2i [TB-[T3| 
and triple well [l3l [l^ potentials yielding a wide variety 
of dynamical scenarios ranging from robust population 
transfer to quantum blocking, non-linear oscillations and, 
even, breaking down of the adiabatic dynamics. 

Following a different perspective, there have been sev- 
eral recent proposals to coherently manipulate single 
atoms [H-iOl and BECs [Hill,!!!] in triple- well poten- 
tials by adiabatically following a particular energy eigen- 
state of the system, the so-called spatial dark state. This 
spatial dark state only involves the two ground states of 
the extreme traps in a close analogy to the well known 
quantum optical Stimulated Raman Adiabatic Passage 
(STIRAP) technique [23[. Accordingly, these tools were 
named Three-Level Atom Optics (TLAO) techniques 
[isj . In this paper, following the work by Vitanov and 
Shore [l^l on the adiabatic passage on two-level atoms, 
we extend the TLAO techniques to the two-level matter 
wave case showing that the adiabatic dynamics of a BEC, 
in a double well potential, can be described in terms of a 
dark variable resulting from the combination of the pop- 
ulation imbalance and the spatial atomic coherence. This 
dark variable can be tailored varying in time the tunnel- 
ing rate, the energy bias, and the non-linearity with the 
goal of (i) adiabatically splitting of the BEC, (ii) achiev- 
ing complete BEC transfer from one well to the other; 
and (iii) preparing of an adiabatic self-trapping state. 

The paper is organized as follows. In Section II we 
present the physical system consisting of a BEC in a 
double-well potential. We assume the two-level approx- 
imation and describe the BEC dynamics in terms of a 
dark variable. The conditions for the adiabatic control 
of tunneling by means of this dark variable are discussed 
in Section III from a non-linear dynamics perspective. 
In Section IV we present detailed numerical simulations 
on the adiabatic splitting, transport, and trapping of a 
BEC. Section V summarizes the main conclusions of the 
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paper and briefly discusses the validity of the two-mode 
approximation. It also presents some possible extensions 
of the present work such as its formulation in second 
quantization. 



II. MODEL 

We consider a BEC trapped in a double well potential 
(Fig. 1(a)), whose dynamics at zero temperature is de- 
scribed by the time-dependent Gross-Pitaevskii equation 
(GPE) {h = 1): 



' dt 
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where V{r^t) is the external trapping potential and the 
non-linearity is given by g = AN'Kas/m, with TV the to- 
tal atom number, the s-wave scattering length, and 
m the atomic mass. Assuming the two-level approxima- 
tion [2^ , the wave function or classical order parame- 
ter of the BEC under study can be written as ip^fyt) ~ 

satisfying |ci|^ -I- \cr\^ — 1. 'pL.Ri'r) accounts for the 
ground state of the corresponding isolated trap. The 
amplitudes cl and cn obey the non-linear two-mode dy- 
namical equations given by; 
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where Ur^l are the atomic self-interaction energies, f2 is 
the tunneling rate between the two wells, and cr^l are 
the on-site energies. In terms of the 4>l,r{^ overlaps, the 
former parameters read: 
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with / dr[4>*(j3j\ = 5ij j = L,R). 

Notice that the two-mode approximation assumes that 
the parameters of the problem (tunneling rate, non-linear 
interaction, and energy bias) can be varied indepen- 
dently. In general, however, the temporal modification 
of any of these parameters will result in the modifica- 
tion of the spatial mode functions affecting the rest of 
the parameters. From the experimental point of view, 
a double-well optical potential can be created with well 
separations on the range of few microns by means of two 
focused laser beams [2g , or by the superposition of a 3D 
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FIG. 1. (a) Two-level model for the BEC in a double-well 
potential being il the tunneling rate and e the energy bias, 
(b) Three-level correspondence of (a) for the density matrix 
variables and coupling strengths, w is the population dif- 
ference, u/2 (v/2) is the real (imaginary) part of the spa- 
tial coherence, and U is the non-linear self-interaction energy. 
(ft=l) 



crossed beam dipole trap with a ID optical lattice Q . In 
the later, the intensity of the standing wave field that cre- 
ates the optical lattice and its displacement with respect 
to the center of the dipole trap could be used to manipu- 
late both the tunneling and the energy bias. In addition, 
the scattering length could be controlled playing around 
of a Feshbach resonance [27| . 

Within the density matrix formalism, assuming Ur ^ 
Ul{= U), and taking e = eR — e^as the energy bias 
between the two wells, the coherent dynamics of the BEC 
in the two- well potential can be written as follows [l^ : 
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where we have introduced the real-valued variables u 
2 Re {aLR}, v = 2 Im {ctl^r}, w = grr - cjll being ct;; 



(7) 



CiC* and 



CiC* with i,j 



L, R the corresponding 



trap populations and spatial coherence (see Fig. 1(b)). 
Note that the conservation of the norm implies vo^ +v? + 
= 1. 

Inspired by the work of Vitanov and Shore [l^l on STI- 
RAP in a system of two internal atomic levels, we extend 
their proposal to the external degrees of freedom of mat- 
ter waves. It is easy to verify that Eqs. ([7]) are equivalent 
to the discrete non linear Schrodinger equation (DNLSE) 
for a BEC in a triple well potential (see, for instance, 
Ref. under the following identifications: u <^ Cr, 

V ^ —iCm, w <^ Cl, e + Uuj 4^ ^mr and 2Q f^LM, 
where Ci with i = L, M, R are the probability amplitudes 
to find the BEC in the left, middle, or right trap, respec- 
tively. fl^M and ilj\ffl' denote, respectively, the tunnel- 
ing interaction between the left-middle and middle-right 
traps. Note however that for the two level system the 
variables u, v and w are real quantities, while for a three 
level system the variables Ci are, in general, complex 
numbers. 

In the TLAO techniques within the three-level ap- 
proximation being the left and right levels resonant, 
the transfer process is based on adiabatically follow- 
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ing one of the three energy eigenstates of the system 
\D{Q)) = cosO|L) — sin0|i?), with the mixing an- 
gle defined as 6(t) = tan^^ [^LM{t)/^MR{t)], and \L) 
and \R) being the ground states of the left and right 
wells. In the two-level system under investigation, the 
analogue state will be given by a combination of the 
population difference and the coherence d{9) = cos 6* x 
w — sinO X u, with the mixing angle given by 6{t) — 
tan-i [2n{t)/ie{t) + U{t)w{t))]. This analogy opens the 
possibility to extend the TLAO techniques to two-level 
systems, by appropriately engineering the time depen- 
dence of the tunneling rate ^{t), the energy bias e{t), 
and/or the non-linear interaction U{t). Note that the 
non-linear interaction parameter can be modified in time 
by the temporal variation of the scattering length Os us- 
ing either magnetic [s^i or optical [Slj Feshbach reso- 
nances or varying the trap frequency, leading to a mod- 
ification of the BEC spatial profile according to Eq. (5). 
For a review on the manipulation of Feshbach resonances 
see fal. 
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with cither the energy bias e{t), or the BEC non-linearity 
preceding the tunneling interaction fl{t), i.e., the 
counter intuitive sequences At^ > or Atjj > will be 
assumed, n^^j/ = 0, ±1 is a switch that takes the value 
for the two-level matter wave analogue of the STIRAP 
sequence and +1 (—1) for the symmetric (antisymmet- 
ric) double-STIRAP sequences. Adiabaticity means that, 
at any time, \0{t)\ should be much smaller than the en- 
ergy separation between the selected cigcnstatc and the 
one energetically closest. For a weak non- linear interac- 
tion, \U{t)/il{t)\ ^ 1, this adiabaticity condition reads 
\0(t)\ < ^4r2(i)2 + {e{t) + U{t)ojy. It is worth to re- 
mark that for [/ 9{t) is not a parameter of the system 
but a dynamical variable since it contains the population 
difference w in its definition. 



III. ADIABATIC CONTROL OF TUNNELING 



In this section we study different scenarios for the co- 
herent control of the external degrees of freedom of a 
BEC by adiabatically following the dark variable d{9). In 
particular, we show how to adiabatically split, transport, 
and inhibit tunneling of a BEC via the matter wave ana- 
logues of both STIRAP and double-STIRAP techniques 
in two-level systems. 

Let us assume that the BEC is initially prepared in 
the left trap with n{t = — oo) = 0. If so, note then that 
9 = meaning w = —1 and u = v = will be the initial 
state. Then, to coherently split the condensate, 9 should 
vary adiabatically from = to = 7r/2 corresponding 
to M = 1 and w = V = (see Fig. 1(b)) by appropriately 
changing fi, e, and U in time. To transfer the BEC from 
the left to the right trap, the mixing angle should be 
slowly increased up to = tt. For the adiabatic self- 
trapping state case, the evolution will consist in varying 
9 from to e.g., tt/2 and then back to 0. In all cases, in 
order to guarantee that the BEC follows the dark variable 
d{9) during the whole process two conditions must be 
fulfilled: 

Condition 1: Adiabaticity criteria, ^{t), e(i), and U{t) 
should be smoothly varied in time to adiabatically follow 
the dark variable d{9) which, in turn, means that the 
mixing angle 9{t) = iim^'^ 2 [VL{t) / {e{t) + U{t)w)] must 
be slowly changed from 9{t = —oo) = to its expected 
final value. Gaussian profiles for the temporal variations 
of the control parameters are assumed: 



n = rig + r2oe-(*-*"' 
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Condition 2: Avoiding adverse bifurcation points. For 
large enough values of the non-linearity, the interaction 
between the atoms of the BEC results in a non-linear 
temporal coupling producing additional non-linear sta- 
tionary states yielding loop structures and a rich variety 
of level crossing scenarios [l^, [H, [H, H^l • In the mat- 
ter wave STIRAP case for a triple- well potential [31, it 
has been shown that even in the adiabatic limit given 
by \9{t)\ — > 0, the appearance of non-linear stationary 
states breaks down, in some cases, the adiabatic evolu- 
tion. To address this issue in the double-well potential, 
we start first looking for those critical U values giving rise 
to non-linear stationary states by solving the eigenvalue 
equation: 



n 



CL 

cr 



CL 



(11) 



with % given in Eq. (3) and /i being the chemical po- 
tential. After some algebra one obtains the following 
fourth-order eigenvalue equation for ^: 



{-2ii + e + 2Uf {-2n + e + Uf 
+'in^{~2fi + e + 2Uf = 



(12) 



For e = and U^/{2fl)^ < 1_(> 1) this quartic equation 
gives two (four) real roots [l^l ■ 

Additional information on the BEC dynamics can be 
obtained looking for the stationary solutions of the den- 
sity matrix equations (7) and analyzing their stability 
(see Appendix A). The stationary solutions can be writ- 
ten as {u^", w^" = 0, w''**} with u^'^/w^^ = tan 6*^" and 
^yss-j2 _^ (iy'ss)2 _ ^ rpj^^ energy of these stationary so- 
lutions, up to four, is given by Eq. (12). Note that, as 
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the dynamics is conservative, the eigenvalues of the hnear 
stabihty matrix must satisfy E^^j^A^ = for each station- 
ary solution. In fact, a Linear Stability Analysis (LSA) 
around any of these solutions yields one eigenvalue equal 
to zero together with either two pure imaginary eigen- 
values corresponding to a stable fixed point (or center) 
or two real eigenvalues accounting for an unstable saddle 
point. In the limit of a large non-linear interaction, i.e., 
for U"^ /[2nY > 1, the system presents four stationary so- 
lutions, one of them being an unstable saddle point while 
the other three are stable elliptical ones. 

If the number of stationary solutions remains constant 
during the whole adiabatic dynamics, the system will suc- 
cessfully follow the selected energy eigenstate. However, 
depending on the interplay between the non-linear in- 
teraction and the tunneling rate, during the dynamics 
the number of stationary solutions will change through 
a bifurcation point from four to two or vice versa. If 
the selected stationary solution is not involved in the bi- 
furcation, the system will again adiabatically follow the 
corresponding eigenstate. On the opposite case, whether 
the adiabatic dynamics will be affected will depend on 
the particular bifurcation scenario: 



Case 1. From 4 to 2 stationary solutions through a back- 
ward pitchfork bifurcation point in which the stationary 
solution to be followed merges two more solutions (one 
unstable) and yields one stable solution. In this case, the 
system will evolve adiabatically. In what follows, this 
case will be termed PB42 case, even if the pitchfork bi- 
furcation is perturbed and becomes imperfect. 
Case 2. From 2 to 4 solutions through a forward pitchfork 
bifurcation point in which the stationary solution to be 
followed becomes unstable and yields two stable solutions 
corresponding, in the limit of vanishing tunneling, to w = 
±1. In this case, the system after reaching the bifurcation 
point will split into a combination of the two stationary 
solutions. In this case, named PB24 in what follows, the 
adiabatic dynamics will, in general, break down. 

Case 3. From 4 to 2 solutions through a saddle-node 
bifurcation in which the stationary solution to be fol- 
lowed is annihilated by an unstable solution. The adia- 
batic dynamics breaks down even in the adiabatic limit. 
After reaching the bifurcation the system will split into 
a combination of two non-degenerate energy cigenstatcs 
and, therefore, will oscillate at the corresponding energy 
splitting. Case SNB42 in what follows. 



IV. NUMERICAL SIMULATIONS 

This section is devoted to the numerical simulations of 
the EEC adiabatic splitting, transport, and self-trapping 
by means of the temporal variation of the energy bias. 
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FIG. 2. (Color online) Adiabatic splitting of a BEC. (a) Pop- 
ulation difference w{tout) at the end of the splitting process 
in the parameter plane modulation amplitude of the energy 
bias 60 versus the non- linear interaction parameter Ug. Initial 
conditions: w{ti„) — —1 and u{tin) = v{tin) — 0. Parameter 
setting: a, = an = 212.8^^ S At, = SOOfio \ tn = 1500^0 ^ 
and Q,g = eg = — Uo = nu = 0- (b) As in (a) for the 
fixed value eo ~ l-5Qo- The temporal variation of the en- 
ergy bias and the tunneling rate is shown in the inset. From 
a non-linear dynamics perspective, we have classified the re- 
sults into five difi^erent regions, from a to e, depending on the 
number and type of bifurcations that the system suffers dur- 
ing its dynamics (see text). The detailed temporal dynamics 
corresponding to cases (i) and (ii) are plotted in Figs. 3 and 
4, respectively. The adiabatic splitting succeeds in region c 
corresponding to Ug/Q,o — (—0.12, 1.5). 



the non-linearity, and the tunneling rate. We will show, 
for the most relevant cases, the eigenvalues and station- 
ary states predicted by Eqs. (12) and (7) as well as their 
linear stability. The main goal of these numerical simu- 
lations will consist in illustrating the different dynamical 
scenarios to adiabatically control the BEC while look- 
ing for those parameter values that prevent the system 
to reach the two unwanted bifurcation cases PB24 and 
SNB42 previously described. 
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FIG. 3. (Color online) Adiabatic splitting of a BEC. Ug = 
O.SSlo and the rest of parameters as in Fig. 2. As a function 
of time: (a) non-linear interaction, energy bias and tunneling 
rate profiles; (b) energy eigenvalues; (c) stationary solutions 
for the population difference w""; and (d) population differ- 
ence w{t) after numerical integration of Eqs. (7). In (b) and 
(c), the dashed curve accounts for the unstable solution, while 
the short arrows indicate the adiabatic solution (in red) to be 
followed by the system. PB42 and PB24 denote the two pitch- 
fork bifurcations that the system suffers at Qot — 777 and 
Qot = 1750, respectively. Note that the pitchfork bifurcation 
PB24 yields one unstable solution (dashed curve) plus two 
energetically degenerated stable solutions (solid curve) corre- 
sponding, in the limit of vanishing tunneling, to lu = ±1. 



A. Adiabatic Splitting of a BEC 

We start first by fixing the non-linear interaction con- 
trol parameter U to a constant value Ug while, counter in- 
tuitively, time-varying the energy bias and the tunneling 
rate (see the inset of Fig. 2(b))) with the goal of achieving 
equal population in the two wells, i.e., w(tout) = start- 
ing with the BEC being in the left trap, i.e., w{tin) = —1, 
or, equivalently, 9[tout) = 7r/2 from 0{tin) = 0. Fig. 2(a) 
shows w{tout), after the numerical integration of Eqs. (7) 
with w{tin) = — 1 and the rest of parameters values given 
in the figure caption. Clearly, there is a large set of 



FIG. 4. (Color online) As in Fig. 2 for Ug = l.rQo- SNB42 
and PB24 account for a saddle-node and a pitchfork bifurca- 
tion at Hot — 1374 and Qot — 1588, respectively. 



parameters where the adiabatic splitting process takes 
place with a high fidelity (see the plateau in Fig. 2(a)). 
Fig. 2(b) shows the plateau w{tout) = for Ug/ilo = 
(—0.12, 1.5) and eo = 1.5r2o- To give more insight into 
the dynamics of the system, we have divided the non- 
linear interaction range studied in five regions, from a 
to e, as shown at the top of Fig. 2(b). In the following 
we will discuss in detail the dynamics for each of these 
regions. Note that in all cases, the dynamics starts and 
ends with fl(tin) = ^{tout) = which, for U ^ 0, implies 
that for t = tin,out the system presents four fixed points 
or stationary solutions. For intermediate times and large 
enough tunneling rates, there are only two stationary so- 
lutions. 

In Fig. 3 we show the temporal dynamics correspond- 
ing to the parameter values of case (i) in region (c) of 
Fig. 2(b) with Ug = O.bflo (see Fig. 3(a)). The adia- 
batic energy eigenvalues and stationary solutions w'"' are 
shown in Figs. 3(b) and (c), respectively, with the dashed 
curve accounting for the unstable ones. Short arrows in 
Figs. 3(b) and 3(c) indicate the stationary solution to be 
adiabatically followed. The dynamical variable w(t), af- 
ter integration of Eqs. (7), is plotted in Fig. 3(d). For 
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this set of parameter values the adiabaticity condition is 
fulfilled and the selected eigenstate is involved in a sin- 
gle bifurcation, at ^ot = 777, corresponding to the PB42 
case previously described. After the bifurcation point, 
the system follows the eigenstate that, at the end of the 
process, yields w = and fj. = 0.25Qo. This non-linear 
dynamical scenario holds for the whole plateau, region c 
in Fig. 2(b), where the adiabatic splitting succeeds. 

For the parameter range Ug/ilo = (1.5,2.07) corre- 
sponding to region d in Fig. 2(b), the final state of the 
BEC strongly depends on the parameter values. For a 
slight modification of Ug, the final state alternates be- 
tween w ^ +1 and w ^ — 1 accounting for the BEC 
being located at the right or left trap, respectively. To 
understand the origin of this behavior, we plot in Fig. 4 
the detailed dynamics for Ug = I.TSIq corresponding to 
case (ii) in Fig. 2(b). At ftot = 1374, the solution that 
the system is adiabatically following is annihilated with 
an unstable one through a saddle- node bifurcation, i.e., 
case SNB42 previously described. At this point, the 
system splits into two non-degenerate components (be- 
ing the largest one the energetically closest stable solu- 
tion) and starts to oscillate at the frequency correspond- 
ing to the energy separation of the two solutions. At 
Ugt = 1588, the system reaches a bifurcation point of 
the PB24 type. This bifurcation yields two energetically 
degenerated stable solutions corresponding to the two fi- 
nal states w — ±1. Thus, at ^Igt = 1588, the largest 
component of the system chooses one out of the two new 
stable solutions. We have numerically verified that the 
selected solution strongly depends on the previous oscil- 
latory dynamics. 

For \Ug\ > 2.07f2o corresponding to regions a and e 
of Fig. 2(b), there are four stationary solutions during 
the whole dynamics that do not present any bifurca- 
tion. In this case, the non-linear coupling forces the adi- 
abatic evolution from 9 {tin) = to 6 {tout) ~ instead of 
evolving to the desired 6 {tout) = 7r/2. Then, the initial 
{u,v,w) = (0,0,-1) and the final state of the system 
coincide. 

Finally, for -2.070o < Ug < -0.1217o, region b in 
Fig. 2, the energy eigenstate to be followed reaches, 
during the ramping down of the tunneling interaction 
51 (t), a pitchfork bifurcation of the PB24 type. After 
this bifurcation point, the two stable stationary solutions 
are energetically degenerated and, therefore, the system 
splits into a non-oscillatory combination of them that 
for Ug = -O.nflo {Ug = -2.0717o) yields w{tout) = 
{w{tout) = -!)■ 

Up to now, we have discussed the possibility to adia- 
batically split a BEC by fixing U and time varying e{t) 
and n(i). Note, however, that it is also possible to adi- 
abatically split the BEC by appropriately time varying 
the non-linear interaction U {t) and the tunneling rate 
n{t) with a fixed energy bias. Thus, for Ug = —OAflo, 
Uq ~ 1.551oi e = and the rest of the parameters given in 
the figure caption. Fig. 5 shows w{tout) as a function of 
the time delay between the modulation of the tunneling 
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5. Adiabatic splitting of a BEC. Population difference 
out) at the end of the BEC splitting process as a function of 
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the temporal delay Atu- The temporal variation of the non- 
linear interaction and the tunneling rate is shown in the two 
upper figures for Atu = —1250,^^ (left) and Atu ~ 425Q. 
(right). Initial conditions: w{tin) = - 
0. Parameter setting: Ug — —OAQq, 
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Uo — 1.50,0, ou — on = 
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adiabatic splitting succeeds in the parameter region f^o At[/ 
(250,600). 



rate and the non-linear interaction. We again obtain a 
plateau of parameter values, ^o/S.tu — (250,600), where 
the adiabatic splitting, corresponding to w{tout) ~ 0, be- 
comes successful. A detailed analysis of the temporal dy- 
namics reveals that, as in Fig. 2, the plateau corresponds 
to the case where the selected solution reaches during its 
dynamics a single bifurcation point of the PB42 type. 



B. Adiabatic Transport 

To adiabatically transfer the BEC from the left to the 
right trap, the mixing angle must be slowly varied from 
0{tin) = up to 9{tout) = TT. Taking rij = — 1 and the rest 
of parameters as in Fig. 2, Fig. 6 shows the population 
difference w(to«t) at the end of the antisymmetric double 
STIRAP process (shown in the inset) as a function of 
the non-linear interaction Ug. Clearly, there is a wide 
parameter setting, the plateau Ug/ilo — (—2, 1.4), where 
the BEC transfer process takes places with a high fidelity. 
In the parameter region Ug/fto = (—2,0), the system 
starts at w{tin) = — 1 following an energy eigenstate of 
the system that does not involve any bifurcation point. 
For Ug/^lo = (0,1.4) the system crosses first a PB42 
bifurcation to later on reaching a PB24 bifurcation point 
that, for these parameters, yields eventually w{tout) ~ 1- 
In regions Ug/Vto = (-4,-2) and Ug/^a = (1-4,4) the 
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FIG. 6. Adiabatic transport of a BEG. Population differ- 
ence witout) at the end of the BEG transport process as 
a function of the non-linear interaction parameter Ug. The 
temporal variation of the energy bias and the tunneling rate 
is shown in the inset. Initial conditions: w{tin) — —1 
and u{tin) = v{tin) = 0. Parameter setting: eo = 1.5Qo, 
(7, = an = 212.8^0 \ AU = SOO^o \ tn = ISOOt^o \ n, = -1 
and Uo = tg = nu = 0. The adiabatic transport succeeds in 
the parameter region Ug/^lo = (^2, 1.4). 
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FIG. 7. Adiabatic self-trapping of a BEG. Population differ- 
ence at (a) the intermediate time t = tn and (b) at the end of 
the BEG transport process, as a function of the energy bias 
eg. The temporal variation of the non- linear interaction and 
the tunneling rate is shown in the inset for Uo — l.Silo. Initial 
conditions: w{tin) — —1 and u{ti„) = v{ti„) = 0. Parameter 
setting: Ug = -0.4^0, au = an = 212.80(7 \ Atu = 325^1^^, 
tsi = 1500r2,7^, nu = 1 and eo — Uo —12^=0. The adiabatic 
self-trapping state does not succeed in the parameter region 
eg/no = (0,0.8). 



C. Adiabatic Self Trapping 

It is also possible to completely inhibit the BEC trans- 
port between the two wells by means of the matter wave 
analogue to the double STIRAP process (see inset of 
Fig. 7). Starting with the BEC located at the left trap, 
w(tin) = —1, we show in Fig. 7 the population differ- 
ence at the time when the tunneling rate is maximum, 
w{tn), and at the end of the double-STIRAP sequence. 



w{to 



as a function of the energy bias €g. The adia- 



non-linearity yields 9{to^^t) 
final state to be the same. 



forcing the initial and 



batic self-trapping process succeeds for a wide domain of 
parameter values, except for the region eg/51o = (0, 0.8). 
In this last region, the solution to be followed reaches 
during its dynamics a saddle-point bifurcation while, for 
the rest of values of the energy bias, the selected solution 
does not involve any bifurcation point. 



V. CONCLUSIONS 

By means of a dark variable that results from the 
combination of the population imbalance and the spa- 
tial atomic coherence, we have investigated in detail the 
adiabatic dynamics of a BEC in a double well potential 
in the framework of the two-level approximation. We 
have shown that it is possible to robustly split, transport 
or trap a BEC by appropriate temporal variation of ei- 
ther the energy bias or the non-linear interaction together 
with the tunneling rate. All these proposals have been 
studied from a non-linear dynamics perspective deriving 
the stationary solutions of the system, evaluating their 
linear stability, and discussing the bifurcation scenarios. 
For the eventual implementation of the techniques here 
discussed, we want to highlight the recent work by M. 
Bauer et al. p7| reporting an accurate control on mag- 
netic Feshbach resonances by means of laser light. 

In order to check the validity of the two-mode ap- 
proximation, we have numerically integrated the one- 
dimensional Gross-Pitaevskii (CP) equation for the non- 
linear adiabatic splitting of a BEC in a double-well po- 
tential, obtaining the characteristic 'plateau' of the STI- 
RAP protocol shown in Fig. 2(b). However, a detailed 
numerical investigation of the GP equation still remains 
to be performed to validate that the matter wave non- 
linear STIRAP techniques here derived in the two-mode 
approximation could be used for matter wave interfer- 
ometry or coherent transport of a BEC. In performing 
this analysis, the global landscape provided by the re- 
sults of the two-mode approximation should be a useful 
roadmap. 

Although it is out of the scope of the present paper, it 
would be very interesting also to extend the present non- 
linear matter wave STIRAP techniques to the second- 
quantization formalism (35j . Within this formalism, we 
could elucidate whether the adiabatic splitting of the 
BEC results in a macroscopic superposition where the 
entire condensate localizes in one of the wells, i.e. in a 
NOON state [S^, or in a perfect 50% spatial splitting of 
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all the atoms. Note that in both cases the final popula- 
tion difference reads w(tout)=(^- 



lutions in the form of: 
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Appendix A: LINEAR STABILITY ANALYSIS 

The stationary solutions of the density matrix 
equations (7) take the form w'*'' = 0, w*''} with 

u''/w'' = 2n/{e + Uw"') and {u'')^ + {w"")^ 



1. Let 



us consider now a perturbation around any of these so- 



u{t) = u'" + 5ue^* 

v{t) = (5we^* 

w{t) = w""" + Swe^* 



(Al) 
(A2) 
(A3) 



where the real part of A determines the stability of the 
solution. The stability is governed by the equation (M — 
XI)6s = where 5s = {6u, Sv, 5w)^ and M is the linear 
stability matrix defined as: 



M = (e 










2n 











2f7 



(A4) 



Solving the secular equation det(A/ — A/) = 0, one ob- 
tains the following three eigenvalues: 



Ao = 



(A5) 



A± = ±yj2Vt{Uu'''' - 2Vt) - (e + C/u;-^'')^ (A6) 

As expected from a conservative system, S^^j^A; = 
for every stationary solution. One eigenvalue is equal 
to zero while the other two are either pure imaginary 
corresponding to a stable fixed point (or center), or pure 
real accounting for an unstable one (or saddle point). 
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